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ABSTRACT 


Dean vortices in curved channels with 40 to 1 aspect ratios were 
measured and studied over the range of Dean numbers from 50 to 450. At 
low Dean numbers (<50) the flow is fully laminar. At higher Dean numbers, 
the development of vortex pairs as the primary instability was a function 
of Dean number and the angle of curvature. Higher Dean numbers required 
less curvature to produce the vortices. At Dean numbers from 75 to 200, 
secondary instabilities developed in the form of vortex pair undulations, 
vortex pair twisting, and in the form of events where vortex pairs appear 
and disappear. These secondary instabilities are also a function of Dean 
number and curvature. Twisting leads to increases in longitudinal 
fluctuating intensities, particularly in the upwash region from the 
concave wall, which are especially significant at Dean numbers above 150. 
Fluctuation increases eventually lead to fully turbulent conditions and 
depend on location in the spanwise/radial plane relative to a vortex pair 
structure, where upwash regions from the concave wall are the most 
unstable. Fully turbulent flow develops at Dean numbers greater than 


about 400. 
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i INTRODUCTION 


The influences of curvature on flow behavior are important to a 
variety of machines including pump impellers, radial compressors, nozzles 
and diffusers, axial flow compressors, and axial flow turbines. To 
develop predictions of these flows, bench mark experimental results are 
required. The present study is part of an experimental effort to achieve 
this aim. In it, flow instabilities in a constant cross section channel 
with constant curvature are examined. These results not only enhance 
understanding of flows in more complex geometries, but also provide a 
building block for the development of improved prediction schemes and more 


efficient machine designs. 


A. BACKGROUND 

Until recently, only two laminar flows were known to occur prior to 
the onset of turbulence in curved channels. At low enough Re, the most 
basic laminar flow occurs. This is curved channel Poiseuille flow (CCPF). 
At higher Re, secondary flow develops due to centrifugal instability of 
CCPF. This second laminar flow consists of an array of streamwise- 
oriented, counter-rotating pairs of vortices and is called Dean vortex 
flow. 

If the spanwise spacing of vortex pairs is A, then a non-dimensional 
spanwise wavenumber of the vortex pairs is a=2(d/2)m/A. For given a and 
radius ratio p, classical linear stability analysis predicts the neutrally 
stable De which sets the lower limit where such vortices may be present. 
The minimum of the neutral stability curve for the primary instability 
then occurs at De, and a. For n=0.979, 16.236 

Recent numerical work by Finlay et al. (Ref. 1) predicts these Dean 
vortices become unstable at higher Re to two types of travelling waves, 


resulting in two kinds of wavy vortex flows, referred to as undulating 


Dean vortex flow and twisting Dean vortex flow. The first mention of 
experimentally observed waviness is given by Kelleher et al. {Ref. 2], and 
the first photographic evidence of vortex motions believed to be twisting 
is presented by Ligrani and Niver (Ref. 3). Additional experimental 
evidence of wavy vortex flow in a rotating, curved channel is given by 
Matsson and Alfredsson [Ref. 4]. Finlay and Nandakumar (Ref. 5] provide 
experimental evidence of undulating Dean vortex flow and twisting Dean 
vortex flow which are compared to wavy vortex results obtained by direct 
simulation of the time-independent  three-dimensional Navier-Stokes 
equations. The experimental results were obtained using the same large 
aspect ratio curved channel described by Ligrani and Niver (Ref. 3). 

Another category of motions which has been observed is related to 
appearances and disappearances of Dean vortex pairs across the span of the 
channel. Evidence of such behavior for channel flow was first provided by 
Ligrani and Niver (Ref. 3] in the form of photographic sequences showing 
small secondary vortex pairs observed temporarily in spanwise/radial 
planes for De ranging from 75 to 220. Additional study of different 
phenomena that result in these appearances and disappearances is presented 


by Ligrani et al. (Ref. 6]. 


B. OBJECTIVES 

The present study focuses on the behavior of Dean vortices An 
rectangular curved channels with 40 to 1 aspect ratios. Methods are 
employed to study flows at De from 49.2 to 450. This range of De 
represents flow regimes from laminar CCPF to fully turbulent flow. Of 
particular interest are undulation, twisting, appearances of vortex pairs, 
disappearances of vortex pairs, and their connections to the development 


of fully turbulent flow. 


C. ORGANIZATION 
Chapter II describes the construction and operation of the curved 


channels used for the present study. Included is information on the 


determination of bulk flow rates in the channels. Chapter III describes 
equipment used for measurement of pressure and velocity surveys plus 
techniques used to obtain real time flow visualization results. Chapter 
IV presents experimental results. Chapter V presents the summary and 
conclusions. Appendix A contains uncertainty analysis results. Appendix 
B describes the software employed. Appendix C contains sequential Dynamic 
Signal Analyzer steps to obtain power spectra plots. Appendix D contains 
all tables. Appendix E contains photographs, schematics and result 


figures. Appendix F contains individual power spectra plots. 


LI. EXPERIMENTAL FACILITIES 


A. CURVED CHANNEL DESCRIPTIONS 

Two nearly identical curved channels, a transparent channel and a heat 
transfer channel, were utilized. A photograph of the transparent channel, 
a photograph of the blower assembly, and a schematic of the channel are 
shown in Figures 1, 2 and 3, respectively. The channel is described in 
detail by Ligrani and Niver [Ref. 3]. The facility is an open-circuit 
suction facility designed for low-speed transition studies. At the inlet, 
an aluminum honeycomb and three screens reduce spatial non-uniformities in 
the flow. These are followed by a two-dimensional nozzle with a 20 to l 
contraction ratio which accelerates the flow into a 2.44 m long straight 
duct. The straight section provides fully developed laminar flow at the 
inlet of the curved section which then follows, for Dean numbers up to 640 
(Ref. 3]. The radii of the concentric convex and concave walls in the 
curved section are 596.9 mm and 609.6 mm, respectively. The 180 degree 
curved section is then followed by another 2.44 m long straight duct, and 
then a honeycomb, several screens, a diffuser and an outlet plenum. From 
the inlet of the first straight section to the outlet of the second 
straight section, the interior dimensions of the facility are 1.27 cm for 
the height and 50.80 cm for the width, providing an aspect ratio of 40 to 
1 د‎ Tolerances for the height and width are 0.015 cm and 0.050 cm, 
respectively, maintained using longitudinal and cross-beam supports along 
the length of the test section. These supports are connected to two 
continuous, seamless sheets of polycarbonate (commercial name is Lexan) 
which form the two primary interior walls adjacent to the flow stream. 
The outlet plenum (labeled #1 in Figure 3) is connected to a second plenum 
(labeled #2 in Figure 3) using 50.8 mm diameter piping. Plenum #2 is 
maintained at low pressure by a 1/3 horsepower (250 W) blower manufactured 


by ICG Industries, Inc. A brass gate valve is used to adjust the flow 


rate, which is determined from the pressure drop across an orifice plate 
located in the pipe between the two plenum chambers. 

In the transparent channel a slot, 3.2 mm wide and 76.2 mm long is 
located in the convex wall. It is aligned in the spanwise direction 50.8 
mm off the centerline at a location 120 degrees from the start of 
curvature. A support block of polycarbonate (152.4 mm by 63.5 mm by 22.2 
mm) is joined using epoxy to the wall surrounding the slot to maintain the 
convex wall's dimensional integrity, and to minimize the deflection of the 
curved wall in the vicinity of the slot. Foam is used to line the slot to 
allow probe insertion with no air leakage to the channel interior. This 
foam is carefully positioned to avoid any bumps or discontinuities in the 
flow. Additional details are given by Ligrani, et al. [Ref. 7] 

In the second curved channel, a slot identical to the one previously 
described is located straddling the centerline, aligned in the spanwise 
direction, and positioned 19.05 cm downstream of the convex surface. This 
channel is configured for heat transfer measurements, but the slot 
location is the only significant difference pertaining to this study. 


Additional details are given by Hughes. (Ref. 8] 


B. DEAN NUMBER DETERMINATION 

Dean number, channel Reynolds number, channel bulk velocity and mass 
flow rate are determined by measuring the pressure drop across the 3.81 cm 
(1.5 inch) diameter ASME sharp edged orifice plate. Also measured are 
barometric pressure, and ambient temperature. The expansion coefficient 
Y is then calculated from Holman and Gajda [Ref. 10]. An initial Dean 
number guess is then made by interpolating a table of Dean number versus 
orifice pressure drop compiled by Niver [Ref. 11]. The density of air is 
computed using the ideal gas law. An initial mass flow rate is 
subsequently calculated based on the initial Dean number. The Reynolds 


number through the pipe is calculated using the equation given by: 


ae 
Re =m 
Fr pA 





Flow coefficient K is interpolated from ASME tables (Ref. 12]. With this, 


mass flow rate is determined using an equation of the form: 


m=KA, Y [20A P, 


Iterations continue until convergence to within 1% of previous mass flow 
rate. Channel bulk velocity, channel Reynolds number and Dean number are 


then calculated using the equations: 
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respectively. A plot of pressure drop versus De is shown in Figure 4 for 
a barometric pressure of 29.92 inches and an ambient temperature of 19 


degrees Celsius. 


III. EXPERIMENTAL PROCEDURES 


A. STREAMWISE, SPANWISE AND RADIAL VELOCITY SURVEYS 
1. Five-hole Pressure Probe Description 

A photograph and two schematic diagrams of the miniature five-hole 
pressure probe are shown in Figures 5, 6 and 7, respectively. The 
diameter of the probe tip is 1.38 mm. The probe is arranged such that the 
central tube is surrounded by four other tubes. The ends of each of the 
surrounding tubes are beveled at a 45° angle with respect to the mouth of 
the central tube. Each of the five-hole tubes is made from stainless 
steel hypodermic tubing with inner and outer diameters of 0.203 mm and 
0.406 mm, respectively. The probe was constructed by first spot welding 
together these five tubes. The tubes were then bent in a 90° turn with a 
radius of 1.52 mm to avoid any tube crimping. Silver solder was then used 
to join the five tubes together continuously along their length. Finally, 
portions of the probe tip were beveled to appropriate angles using a 
surface grinder. This approach is slightly different than the one 
described by Treaster and Houtz. [Ref. 13] 

Figure 6 shows that the present probe tip extends 5.08 mm (4 tip 
diameters) beyond the bend in the stem to minimize the effects of stem 
flow blockage at the measuring location. Similarly, the tip is 22.9 mm 
(19 tip diameters) from the 9.53 mm diameter mounting cylinder. Within 
this cylinder, each tube is connected to a 1.60 mm diameter tube using a 
small manifold for each tube connection. The mounting cylinder was sized 
for easy handling and mounting. On top of the cylinder is a second set of 
manifolds with five male connectors for attachment to 6.35 mm outer 
diameter/4.67 mm inner diameter Imperial Eastman Polyflo tubing. 

The miniature five-hole probe was designed to minimize the adverse 
effects of flow blockage and spatial resolution when used for measurements 


in the curved channel. When in place for measurements at the centerline 


of the channel, the mounting cylinder is not exposed to the flow, and the 
probe produces 4.8% blockage over a 12.7 mm x 12.7 mm area. A five-hole 
rather than a four-hole design was chosen since it was easier to calibrate 
and because corrections for spatial resolution are more readily applied 
(Ref. 14]. The tip was conical rather than prismatic to minimize near- 
wall flow blockage effects, which are particularly severe for prismatic 
type probes. 
2. Experimental Details 

Figure 8 shows a schematic of the measurement chain used for the 
five-hole pressure probe. Each tube of the five-hole probe is connected 
via flexible plastic tubing (less than 20 cm in length to minimize the 
volume of air in the tubes) to a Validyne model DP103-06 variable 
reluctance differential pressure transducer. A separate transducer is 
used to measure pressure relative to the atmosphere for each port of the 
five-hole probe. Each transducer has a full-scale range of 2.5 mm of 
water differential pressure. Shielded cables connect the transducers to 
separate Celesco CDO1D carrier demodulator units. The carrier demodulator 
units produce DC voltage output proportional to the pressure signals from 
the transducers. Each carrier demodulator output is connected to a 
Hewlett Packard (HP) 3498A extender which is controlled by an HP 3497A 
data acquisition/control unit. The HP 3497A is controlled by an HP Series 
300, Model 236 computer which also processes the data. Voltage samples 
are recorded sequentially from the five channels using an overall sampling 
rate of about 10 Hz. 

3. Probe Calibration 

The probe was calibrated by Fields (Ref.15] using a procedure 
similar to Ligrani, et al. (Ref. 7). During calibration, the probe was 
mounted in a calibration sled and placed in an open-circuit blower tunnel 
which provided uniform flow in the test section. A manually operated 
traverse was used to position the probe through a five yaw angle and five 


pitch angle array. Total pressure was measured with a Kiel probe, and 


tunnel static pressure was measured near the probe position using wall 
Static taps. For each probe position during calibration, yaw, pitch, 
total and total minus static pressure coefficients were calculated. 
4. Curved Channel Procedures 
a. Traversing Mechanism 

When measurements are made in the channel, the probe is 
oriented using a mount which allows initial yaw and pitch angles to be 
adjusted within +0.1° and +0.5°, respectively. The mounting block is 
attached to a travelling block which is part of an automated two- 
dimensional (X-Y) traversing mechanism. Moving parts are driven by 20 
thread per inch (0.8 thread per mm) drive screws which are rotated by 
Superior Electric type M092-FD310 stepping motors. These motors are 
controlled by a MITAS Motion Controller, which directs movement of the 
probe according to instructions from the computer. 

b. Measurements 

The probe is positioned and aligned before measurements are 
made. The probe is first inserted in the slot of the insulated channel at 
a spanwise location 25.4 mm from the centerline so that it just touches 
the straight wall farthest from the slot. It is then visually aligned 
with the flow and then moved to the channel position at the midpoint of 
the spanwise-radial plane using the MITAS controller. Small adjustments 
are made to the yaw and pitch angle by adjusting the probe mount. 
Software program ORIENT, described in Appendix B, is used to check yaw and 
pitch coefficients as this is undertaken until the yaw coefficient is zero 
and the pitch coefficient is sufficiently small. The final orientation of 
the probe gives a yaw coefficient of 0.0+0.01 and a pitch coefficient of - 
0.4. This is equivalent to +3.0° offset. More accurate pitch alignment 
forces the probe against the edge of the slot, hence the pitch offset is 
accounted for during data reduction. No offset 1s applied or needed for 


the yaw angle. 


In a typical spanwise/radial traverse, the probe moves in the 
radial direction from the concave wall to the convex wall as measurements 
are made. To begin measurements at a new spanwise location, the probe 
first moves 2.54 mm in the spanwise direction, then in the reverse span 
direction 1.27 mm, and finally toward the concave wall until it is 1.27 mm 
away. A new radial or spanwise measuring traverse then begins. This 
procedure reduces the adverse effects of backlash and gives repeatable 
data regardless of the direction of spanwise probe motion. 

When measurements are made at a particular location, 50 
samples of each pressure are averaged. This is done by sampling pressure 
from probe ports in sequence so that a long averaging time is provided for 
each channel, thereby minimizing the influence of brief unsteadiness in 
the channel. For each probe location, about 3 min of acquisition and 
processing time is required. A traverse over an entire plane requires 
about 8 h. The entire process is computer controlled and generally 
performed at night when there iS minimal unsteadiness at the channel inlet 
from laboratory disturbances. 

c. Data Reduction 

Measured pressures from each probe port are stored by the 
computer during each traverse. These data are then corrected to account 
for the effect of finite spatial resolution using procedures described by 
Ligrani et al. (Ref. 14]. The procedure corrects pressures to account for 
the fact that they are not measured at the same physical location by 
adjusting pressures from ports 2, 3, 4 and 5 so that they appear to be 
measured at port 1. The correction primarily affects secondary flow 
velocity vectors in the plane normal to the bulk flow direction. With the 
corrections, these secondary vector magnitudes may change as much as 15% 
at locations where velocity gradients are largest. 

Yaw and pitch pressure coefficients are then calculated for 
each location. Using these coefficients and measured pressures, pitch and 


yaw angles are calculated at each location. The total minus static 
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pressure coefficient is then calculated. From this total velocity 


magnitude is calculated using the equation: 


P,-P 
p 


U- AC ots 
The three velocity components are subsequently given by: 

Ug-Ucos(y)cos(p) 

U,=Usin(y) 

and U,=Ucos(y)sin(p) 
where y and p are the yaw and pitch angles, respectively. Mean vorticity 
components are then calculated using finite difference schemes applied to 
the three components of mean velocity. Equations for the vorticity 


components as functions of velocity are given by: 
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Streamwise derivatives are neglected when calculating radial and spanwise 
vorticities. In the final phase of data reduction, corrections are made 
to account for deflection of flow streamlines near the probe tip in the 
presence of transverse gradients of the streamwise velocity. This is 
referred to as the downwash velocity correction [Ref. 14], and is based on 
the idea that the downwash velocity in a particular direction scales with 


the transverse gradient of streamwise velocity in the same direction. 


DM 


B. CONVEX WALL TURBULENCE SURVEYS 
1. 12 Element Hot-film Probe Description 
A photograph of the hot-film array is shown in Figure 9. The 
array has 12 hot-film anemometers mounted equidistant from each other on 
a specially designed aluminum frame. Probe locations are tabulated in 
Table 1, which indicates that probe 1 is positioned 56.95 mm from the 
channel spanwise centerline. Each probe is silver plated platinum with a 
diameter of 0.625 microns and an active length of 0.102 mm. The frame is 
designed to fit snugly into the slot located 120° from the start of 
curvature in the transparent curved channel. With thin shims placed under 
the ends of the frame, the probe is mounted so that the measuring plane of 
the gages is flush with the convex wall. 
2. Experimental Details 
Figure 10 shows a schematic of the measurement chain, including 
the probes, and data acquisition and processing equipment. With the 
present setup, any Six probes can be connected to each of six different 
DANTEC 56C17 Constant Temperature Anemometer bridge units shown in the 
photograph Figure 11. Bridge outputs are sent simultaneously to a HP 
3497A Data  Acquisition/Control Unit and to DANTEC 56N20 Signal 
Conditioning Units. In the discussion which follows, the former signals 
are referred to as raw signals, whereas the latter are referred to as 
conditioned signals. The output of each signal conditioning unit is 
directed to the HP 3497A Data Acquisition/Control Unit which is controlled 
by an HP Series 300, Model 236 Computer which also processes the data. 
With the control unit and computer, voltage samples are recorded 
sequentially from the six channels using an overall sampling rate of about 
5:9 2 
3. Curved Channel Procedures 
a. Measurements 
The data acquisition system is arranged so that the signal 


from the hot-film gage closest to the spanwise centerline is sampled 
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first. Signals from adjacent gages across the span of the array are then 
sampled in sequence. The sequence is repeated until 1000 sampled voltages 
for each gage are obtained. Each bridge uses an overheat ratio setting of 
1.6. Each signal conditioning unit is set with filters to pass signals 
with frequencies from 0.1 Hz to 1.0 kHz. A gain setting of 500 is also 
used. 
b. Data Reduction 

After sampled voltages are stored in the computer, they are 
reduced by a factor of 500 to account for the amplifier gain. Mean values 
of the raw and conditioned voltage signals are then determined. In 
addition, rms voltage for each conditioned channel is also determined. 
Normalized rms voltage for each channel is computed as the ratio of rms 
voltage to the raw mean voltage. 

Covariances are determined for the conditioned voltage signals 


using the equation: 


n 
COV= + y (Eref ,-Eref) (Eobs,-Eobs) 
1-1 


Here, n is the number of samples, ref is the reference channel and obs is 


the observed channel. This is equivalent to: 


e (x) e (x+x?) 


Two different probe configurations were utilized. First, covariances for 
probes 2, 4, 6, 8, 10 and 12 with respect to probe 2 located at z/d=4.73 
were determined for De from 50.25 to 399.8. This sequence is termed ‘in 
to out'.  Covariances for probes 12, 10, 8, 6, 4 and 2 with respect to 
probe 12 located at z/d=7.23 were determined for De=50.25 to 399.8. This 
sequence is termed ‘out to in'. Second, ‘in to out' covariances for 
probes 1 through 6 with respect to probe 1 located at z/d=4.48 were 
determined for De from 49.2 to 400.3. ‘Out to in' covariances for probes 


6, 5, 4, 3, 2 and 1 with respect to probe 6 located at z/d=5.73 were 
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determined for De=49.2 to 400.3. After determining covariances, spatial 
correlations are determined using the equation: 


COV 


COR2 —————————————— 
(rmsEref) (rmsEobs) 


This equation is equivalent to: 


e (x) e (x«xl) 


J/e' (x)? Je (x«x!)? 


Length scale, representing the area under a curve of correlation versus 
z/d is determined by the equation: 
T] 


LS= Y” (COR) (AZ) 
SES 


or alternatively: 
L,= f CORdz 


Here, AZ is the spanwise distance between active probes. The summation 
was used to calculate length scales and is an approximation to the exact 


integral method. 


C. FLOW VISUALIZATION 

In order to obtain information on secondary flow motions resulting 
from the vortices, air in the transparent channel was partially 
contaminated with smoke and then photographed. The device used to 
generate the smoke is described by Ligrani and Niver. [Ref. 3] Smoke 
patterns were illuminated in spanwise/radial planes using spotlights 
directed at spanwise slits through black paper used to line the convex 
surface exterior of the channel. A Sony DXC-M3 video camera with a 
Fujinon-TVZ 1:1.7/10-140 mm DCL-914BY zoom lens was used to record the 
illuminated, time-varying patterns as shown in Figure 12. The camera was 


connected to a Sony VO-6800 portable videocassette recorder which collects 
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images at rates as high as 60 frames per second. Recorded images are then 
played back using a Sony VO-5800 videocassette recorder connected to a 
SONY PVM-1910 Trinitron color monitor. By employing the scan/stop action 
feature of this recorder, photographs of instantaneous images using the 
setup shown in Figure 13 were obtained and then arranged in sequences to 
show time-varying phenomena. New time-varying sequences are presented in 
Chapter IV. This chapter additionally gives results from analysis of 
video tapes previously obtained by Ligrani et al. (Ref. 6]. Data were 
analyzed at De from 60 to 150 at angles from the start of curvature from 
85° to 135°. Of primary interest were frequencies of different appearance 
and disappearance events described by Ligrani et al. (Ref. 6]. 
Tabulations of results are given in Tables 2 and 3. In making these 
tabulations, vortex pairs which moved into or out of the field of view 
during the time sequences were ignored. All frequencies were also 


converted to non-dimensional form using the equation: 


11 
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These are also tabulated in Table 4. 

The two new time-varying sequences presented here were obtained using 
the setup of Figure 13 using a Canon A-1 SLR camera with a 35-70 mm, f3.5- 
22 lens. Settings of 55 mm, f5.6 and speed 1/8 second for Kodak Tri-X 
pan, 400 ASA film were utilized following the recommendations of Longest 
(Ref. 16). 

D. DEAN NUMBER SURVEYS OF POWER SPECTRA AND 
LONGITUDINAL VELOCITY FLUCTUATIONS 
1. Experimental Details 
A DANTEC 55P14 single hot-wire probe is used to obtain surveys of 
power spectra and longitudinal velocity fluctuations. The probe sensor is 
platinum plated tungsten wire with a diameter of 5 microns and a sensor 
length of 1.25 mm. Because support prongs are bent, they are 


approximately parallel to the flow direction at the sensor location, which 
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places the sensor normal to the bulk flow direction. During measurements, 
the probe is rigidly mounted in an X-Y traverse such that the sensor 
itself is oriented in the spanwise direction normal to the flow at a 
streamwise location 120° from the start of curvature. Surveys were 
obtained at each of five locations. These locations are listed in Table 
5 and shown in Figure 16 relative to vortex pair structures. The probe is 
operated at an overheat ratio of 1.8 using a DISA 55M10 Constant 
Temperature Anemometer bridge. The output from the bridge is routed 
through a DANTEC 56N20 Signal Conditioning Unit. Settings on the signal 
conditioner were varied with Dean number range as shown in Table 6. The 
output from the signal conditioner was routed to a DANTEC 56N22 Mean Value 
Unit set at a 10 second integration interval. The output from the signal 
conditioner was additionally sent to a Hewlett Packard 6944A 
Multiprogrammer. This high speed data acquisition system was controlled 
by a Hewlett Packard 300 Series Model 9000 computer which also stored and 
processed the data. Acquisition sampling rate was 2500 Hz if De was less 
than 275 or 5000 Hz if Dean number was higher than 275. 

Another output from the bridge was directed to a second signal 
conditioner operated with amplifier gain of 500, low pass filter of 1.0 
kHz, and a high pass filter of 0.1 Hz. The output from this signal 
conditioner was then sent to a DISA 56N25 RMS Unit set with a 10 second 
integration interval. 

After the DISA 56N20 Signal conditioner, signals were also sent to 
the Hewlett Packard 3562A Dynamic Signal Analyzer. This device is used to 
obtain spectra of the time-varying signals from the hot-wire bridge. 
Figure 14 shows the setup for measurement, data acquisition and 
processing. Figure 15 gives an overview of the internal measurement chain 
of the dynamic signal analyzer. Step by step details are given for setup 
and use of the dynamic signal analyzer for spectral measurements in 


Appendix C. 
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With settings on the signal conditioning units adjusted as shown in Table 
6, frequency spectra were obtained between 0.1 and 25 Hz and between 0.1 
and 10000 Hz. 
2. Probe Calibration and Measurement Procedures 
The probe was calibrated in the test section of the wind tunnel 
located in the Department of Mechanical Engineering laboratory of the 
Naval Postgraduate School. From the calibration, N=0.45, B=3.750916, and 


EOC=2.653925, where the equation given by: 


1/n 
27 2 
سا‎ 
B 


is valid when the sensor is subjected to forced convection. For Re, <0.07 
mixed convection affects sensor response and the relationship between 
velocity and voltage is given by the equation: 

Y=AO+AL*X+AZ*X2+A3*X) 
where  Y-U*, xX=E*-EOC*, AO=-0.68058,  Al1=0.99448, A2=-0.24088 and 
A3=0.024986. 

With the experimental arrangement, longitudinal velocity 
fluctuation magnitudes are determined from digital processing of hot-wire 
bridge voltage signals and from signals read from the DISA 56N25 RMS unit. 
For the former, 20000 voltage samples were obtained. For the latter, the 
normalized streamwise velocity fluctuation magnitudes are calculated using 


the equation from given by: 


E,e' 
spe 


(E, -4) 
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(Ref. 17]. 


17 


IV. EXPERIMENTAL RESULTS 


Surveys of total pressure, mean velocity, velocity perturbation and 
vorticity components are presented as measured using the miniature five- 
hole pressure probe for De from 50.4 to 424.8. Each survey consists of a 
320 point spanwise/radial plane array in the second straight section 
beyond the 180° end of curvature. Thus, the surveys are obtained at a 
streamwise location downstream of the curved portion of the channel at a 
location not subject to the immediate influence of centrifugal 
instabilities. The array covers equispaced points with radial y/d limits 
from 0.1 to 0.8 and spanwise z/d limits from -2.0 to 1.9 where y/d=0.5 and 
z/d=0.0 are the coordinates of the center of the spanwise/radial plane. 

Surveys of spatial correlations versus z/d are presented from six 
point spanwise samples obtained at y/d=1.0 and 120° beyond the start of 
curvature. Spatial correlations provide information regarding 
correlations between different parts of vortex pair structures. Six 
probes were used in two different configurations with probe spacings of 
6.35 mm and 3.15 mn. For the first configuration, every other probe 
element was used, starting with probe 2 for De from 50.25 to 399.8. For 
the second configuration, adjacent probe elements were used, starting with 
probe 1 for De from 49.2 to 400.3. 

Power spectra, streamwise velocity and longitudinal turbulence 
intensity surveys are also presented as measured using the hot-wire probe 
for De from 50 to 450. Each survey is obtained in the spanwise/radial 
plane 120° after the start of curvature. For each De, five probe 
positions were surveyed. The five spanwise/radial plane positions were 
chosen to examine different parts of vortex pair structures. Each probe 
position is listed in Table 5 and shown in Figure 16. 

Frequencies of events where vortex pairs appear and disappear are 
presented for De from 60 to 150. Time lines showing vortex pair 


appearance and disappearance event occurrences are presented for De=100 
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and curvature of 115° on the right side of the channel. Still photography 
timed sequences showing a collapsing occurrence and a splitting occurrence 
are presented also. For each De, these data were obtained from 
Spanwise/radial plane observations which were video tape recorded at 


angles from 85° to 135° from the start of curvature. 


A. SURVEYS OF TIME-AVERAGED QUANTITIES 

Surveys of time-averaged quantities determined from five-hole pressure 
probe measurements are presented in figures 17 through 34. Figures 17 
through 19 present total pressure contours. Figures 20 through 22 present 
streamwise velocity contours. Figures 23 through 25 present streamwise 
velocity perturbations. Figures 26 through 28 present streamwise 
vorticity components. Figures 29 through 31 present radial vorticity 
components. Figures 32 through 34 present spanwise vorticity components. 

l. Total Pressure 

In Figures 17 through 19, surveys of Pamb-Ptotal (in Pascals) 

representing ambient pressure minus measured total pressure are given. 
The higher the index in the flow, the higher the difference in pressure, 
i.e., the higher the vacuum measured. Pressure deficits near the concave 
surface in the form of humps are upwash regions between counter-rotating 
vortex pairs. At De=75 the flow shows evidence of portions of at least 
five deficit regions. At De=125.1 the number of deficits is reduced to 
two. The highest deficits occur at De=175.2. At De>300, the surveys show 
reasonable spanwise uniformity. 

2. Streamwise Velocity 

Surveys of streamwise velocity normalized by bulk velocity are 

displayed in Figures 20 through 22. These show velocity deficits near the 
concave wall which are nearly identical qualitatively to the total 
pressure deficits in Figures 17 through 19. These deficits are less 
intense than ones measured in the curved portion of the channel by Fields 
(Ref. 15] because the mechanism to produce the vortices and accompanying 


secondary flows (centrifugal instabilities) is not present. 
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3. Velocity Perturbation 
In order to show streamwise velocity perturbations, curved channel 
Poiseuille flow is subtracted from streamwise velocity. Such 
perturbations are presented in Figures 23 to 25. At De=100.3 a regular 
perturbation pattern is present across the span of the measurement plane. 
The perturbations become spanwise uniform for Dez300. 
4. Vorticity Components 
Vorticity components in Figures 26 through 34 have positive and 
negative signs. Positive vorticity is contained within the solid contour 
lines while negative vorticity is contained within the dashed contour 
lines. As the plots are arranged so that streamwise flow moves into the 
page, positive streamwise vorticity is clockwise. 
a. Streamwise Vorticity 
Streamwise vorticity contours are shown in Figures 26 through 
28. The streamwise vorticity distributions correlate nicely to total 
pressure and streamwise velocity such that velocity/pressure deficits are 
present between the two vortices within each counter-rotating pair. 
Counter-rotating vortex pairs are clearly evident at the locations of the 
two largest total pressure deficits for De=100.3. At De=175.2 two pairs 
of well defined vortices again correspond to the locations of the velocity 
deficits. At De greater than about 300 streamwise vorticity contours tend 
toward spanwise uniformity. 
b. Radial Vorticity 
Radial vorticity contours are shown in Figures 29 through 31. 
As De increases, the flow becomes spanwise periodic at De=75, until De2300 
when the flow becomes spanwise uniform. 
c. Spanwise Vorticity 
Spanwise vorticity contours are shown in Figures 32 through 


34. These contours are consistent with the total pressure contours and 
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contours of streamwise velocity, where negative vorticity components 
closely match the velocity/pressure deficits. Flow is spanwise uniform at 
De=50.4 and for De2300. 

B. TIME-AVERAGED AND SPATIALLY RESOLVED FLOW 

PROPERTIES NEAR THE CONVEX SURFACE 

Results presented in this section were obtained using an array of hot- 
film gages mounted on the convex surface of the channel (y/d=1.0) 120° 
from the start of curvature. Figure 35 shows the locations of the 12 hot- 
film elements with respect to vortex pairs. Table 1 gives probe locations 
with respect to individual vortices. One expects results from probes near 
Similar vortex characteristics to be correlated with each other. Figures 
36 through 43 give correlations versus z/d referenced to probe 2 for 
probes 2, 4, 6, 8, 10 and 12, for De from 50.25 to 399.8. Figure 44 is 
then a 3-D composite of Figures 36 through 43. Figures 45 through 52 give 
correlations versus z/d referenced to probe 12 for probes 2, 4, 6, 8, 10 
and 12, for De from 50.25 to 399.8. Figure 53 is then a 3-D composite of 
Figures 45 through 52. 

Referring to Figures 44 and 53 shows that, in all cases, correlation 
magnitudes decrease as the distance from the reference probe increases. 
Reference probes are probe 2, located at z/d=4.73, for Figure 44 and probe 
12, located at z/d=7.23 for Figure 53. In Figure 44, sign variations of 
correlations versus z/d are evident at De from 50.25 to 299.6. These 
variations give evidence of regions of positive and negative correlation 
between different parts of vortex pair structures. For De=100, probe 2 is 
located directly over a downwash region. Probe 4, the nearest active 
probe, is located directly over an upwash region and its signal correlates 
negatively with the signal from probe 2. At De2299.6, correlations have 
fully turbulent characteristics, dropping to values which do not vary 
Significantly with z/d. 

Similar qualitative trends are evident in Figure 53. For De=100, 


probe 12 is located over a region of negative vorticity. Probe 10, the 
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nearest active probe, is located over a region of positive vorticity and 
its signal correlates negatively with the signal from probe 12. At 
Dez250.4, correlations have fully turbulent characteristics, dropping to 
values which do not vary significantly with z/d. 

Figure 54 shows the average (over all gages employed) normalized rms 
voltage versus De from 50.25 to 399.8. This curve indicates an increase 
in Vrms/Vmean as De increases from 150 to 200. At De higher than 300, the 
variations with De have lower slopes. 

Figure 55 shows scaled, normalized rms voltage versus z/d for each De 
surveyed. The normalized rms voltage for each channel is scaled by the 
normalized rms voltage measured at De-399.8. Thus the scaled values at 
this De are equal to 1.O. These curves show that spatially varying 
structures occur at De=149.9, 200.7, 250.4, 299.6 and 350.5. 

Length scales shown in Figure 56 represent integrations of the 
correlations with respect to z/d. These evidence spatially varying 
structures between De-100 and De-250. 

The previous results were discussed for an active probe spacing of 
6.35 mm. For a second set of measurements, active probe spacing was 
reduced to 3.18 mm. Results are now presented which were obtained using 
the array of hot-film gages with the smaller probe spacing. Figures 57 
through 64 give correlations versus z/d referenced to probe 1 for probes 
1, 2, 3, 4, 5 and 6, for De from 49.2 to 400.3. Figure 65 is then a 3-D 
composite of Figures 57 through 64. Figures 66 through 73 give 
correlations versus z/d referenced to probe 6 for probes 1, 2, 3, 4, 5 and 
6, for De from 49.2 to 400.3. Figure 74 is then a 3-D composite of 
Figures 66 through 73. 

Referring to Figures 65 and 74 shows that, in all cases, correlation 
magnitudes decrease as the distance from the reference probe increases. 
Reference probes are probe 1, located at z/d-4.48, for Figure 65 and probe 
6, located at z/d-5.73 for Figure 74. In Figure 65, sign variations of 


correlations versus z/d are evident at De from 49.2 to 200.9. These 
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variations give evidence of regions of positive and negative correlation 
between different parts of vortex pair structures. For De=100, probe 1 is 
located directly over a region of positive vorticity. Probe 2, the 
nearest active probe, is located directly over a downwash region and its 
signal correlates negatively with the signal from probe 1. At De2249.6, 
correlations have fully turbulent characteristics, dropping to values 
which do not vary significantly with z/d. 

Similar qualitative trends are evident in Figure 74. For De=100, 
probe 6 is located over a downwash region. Probe 5, the nearest active 
probe, is located over a region of positive vorticity and its signal 
correlates negatively with the signal from probe 6. At De2350.2, 
correlations have fully turbulent characteristics, dropping to values 
which do not vary significantly with z/d. 

Figure 75 shows the average (over all gages employed) normalized rms 
voltage versus De from 49.2 to 400.3. This curve indicates an increase in 
Vrms/Vmean as De increases from 200 to 250. At De higher than 250, the 
variations with De have lower slopes. 

Figure 76 shows scaled, normalized rms voltage versus z/d for each De 
surveyed. The normalized rms voltage for each channel is scaled by the 
normalized rms voltage measured at De=400.3. Thus, as for Figure 55, the 
scaled values at this De are equal to 1.0. These curves show that the 
spatially varying structures occur at De-100.1, 150.7, 200.9, 249.6, 299.2 
and 350.2. 

Length scales shown in Figure 77 represent integrations of 
correlations with respect to z/d. These evidence spatially varying 
structures between De=100 and De=200. 

C. SPECTRAL AND LONGITUDINAL VELOCITY FLUCTUATION 

VARIATIONS WITH DEAN NUMBER 

Figure 16 shows the five hot-wire probe locations where surveys were 
taken in the spanwise/radial plane 120° from the start of curvature. At 


each of these locations, time-averaged streamwise velocity and 
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longitudinal turbulence intensity were measured using a single hot-wire 
probe at De from 50 to 450. The five locations were chosen to obtain 
information on different parts of the vortex pair structure. As shown in 
Figure 16, position 1 is located near the concave wall at the base of an 
upwash region (with respect to the concave wall). Position 2 is located 
in an upwash region at the channel radial centerline. Position 3 is 
located near the concave wall near the base of a downwash region (with 
respect to the concave wall). Position 4 is located near the concave wall 
directly under a vortex. Position 5 is located in a downwash region at 
the channel radial centerline. 

The streamwise mean velocity and longitudinal turbulence intensity 
were each determined two ways. With a digital method, the instantaneous 
streamwise velocity was calculated for each of 20000 samples. The average 
of the 20000 instantaneous velocities is then the mean streamwise 
velocity. The longitudinal velocity fluctuation intensity was then 
determined as the root-mean-square of the instantaneous longitudinal 
velocities. With an analog method, mean velocity and longitudinal 
turbulence intensity were calculated based on mean and rms voltages 
integrated over 10 second time intervals. 

Figures 78 and 79 show plots of the average streamwise velocity using 
digital and analog techniques, respectively, versus De for the five probe 
positions listed in Table 5. The digital and analog streamwise velocities 
coincide nicely. All positions show a velocity increase that is generally 
linear with De. The behavior of the velocity for positions 1, 3 and 4 
varies because spatial locations of upwash regions vary with De. The 
proximity of positions 1, 3 and 4 to the concave wall causes associated 
velocity magnitudes to be generally lower than for positions 2 and 5. 

Figures 80 and 81 show longitudinal turbulence intensities determined 
uSing digital and analog techniques, respectively, versus De for the five 
probe positions. The digital and analog data coincide with each other 


reasonably well. For ۲۵۶160, the intensities are very small for all five 


24 


probe positions . At De higher than 160, intensities generally increase 
with De approximately linearly. Positions 2 and 5 have lower intensities 
because they are farther away from the effects of the concave surface than 
positions 1, 3 and 4. 

Normalized longitudinal turbulence intensities are shown in Figures 82 
and 83. Digital and analog data coincide well at most all De, except for 
values of De=50 when the digitally obtained values are incorrectly 
indicated to be too high. For De higher than 160, longitudinal turbulence 
intensity increases with De. The variations at position 4 are due to 
time-averaged spatial location variance of the upwash regions with De. 
Over this range of De, values for probe positions 2 and 5 are lower than 
others because positions 2 and 5 are farther away from the concave wall. 

Figure 84 gives settings used on the dynamic signal analyzer to obtain 
power spectra over frequencies from O.1 to 25 Hz. Summaries of observed 
local frequency maxima are given in Figures 85 through 89 for De=50 to 
De-200 for each of the five probe positions 1 through 5, respectively. 
Frequencies of 1 Hz to 2 Hz, 5 Hz to 6 Hz, 8 Hz, and 20 Hz to 21 Hz exist 
for all five probe positions. An additional frequency of about 3 Hz 
exists for probe positions 3, 4 and 5. Figures 90 and 91 are two examples 
of power spectra used to obtain such information. Figures 85 through 89 
were compiled from data presented in Figures 92 through 151, which show 
individual frequency spectra for all five probe positions at De from 50 to 
450. 

Figure 152 gives settings used on the dynamic signal analyzer to 
obtain high frequency power spectra over frequencies from 0.1 Hz to 10.0 
kHz. Figure 153 shows the average of five power spectra at probe position 
2 at De=170. Clearly evident is a principal peak between 90 Hz and 100 
Hz. Also apparent are harmonic peaks at about 180 Hz, 280 Hz and 400 Hz. 
Such variations are believed to be due to twisting (Ref. 18]. Figure 154 
shows a plot of principal peak magnitudes versus De for all five probe 


positions at De from 150 to 200. In this figure, peak magnitudes are 
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shown to be dependent on De and probe position, where peak magnitudes 
generally increase with increasing De. The highest peak magnitudes 
observed are for De=180 at probe position 4 and for De=200 at probe 
position 1. Such behavior suggests higher instability near the upwash 
region close to the concave surface. Figure 154 was compiled from data 
shown in Figures 155 through 214 which are individual power spectra for 
all five probe positions at De from 50 to 450. 

D. ANALYSIS OF APPEARANCE AND DISAPPEARANCE PHENOMENA 

FROM FLOW VISUALIZATION RESULTS 

Appearances and disappearances of already developed vortex pairs along 
with accompanying alterations of spanwise vortex pair spacing are 
discussed by Ligrani et al. [Ref. 6] At the time of this writing, there 
is evidence of two types of appearance events and four types of 
disappearance events. Appearance events are categorized as growth or 
splitting. Growth events occur when vortex pairs emerge from the flow 
region near the concave wall. Splitting events occur when one vortex pair 
is rearranged into two vortex pairs. Splitting is also observed after a 
vortex pair decreases in size with simultaneous emergence of two new 
vortex pairs nearby from flow near the concave surface. 

Disappearance events are categorized as engulfment, cancellation, 
merging and collapse.  Engulfment refers to the absorbtion of a vortex 
pair by another vortex pair where contact occurs between vortices with the 
same sign of vorticity. Cancellation occurs when vortices of the opposite 
Sign collide. Typically, these vortex pairs rock slightly toward each 
other as they move together. With merging, a smaller vortex pair merges 
into the upwash region of a larger pair. Collapse events are evidenced by 
the disappearance of a pair into the flow region near the concave surface. 

Figure 215 shows a new type of splitting event which was unrecognized 
prior to the present work. Here, the center vortex pair in the field of 
view shrinks in size until it disappears into the flow near the concave 


Surface while two new vortex pairs emerge simultaneously from the base of 
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the shrinking vortex pair. The sequence starts at the top with frame 1. 
Each photo shows the full channel height, 12.7 mm, and a channel width of 
50.8 mm, and gives a spanwise/radial plane view with the concave surface 
at the bottom and streamwise flow into the plane of the page. In frames 
1 through 5 the center vortex pair begins to shrink. In frame 6, dark 
spots begin to form on each side at the base of the shrinking vortex pair. 
Both spots rapidly grow in frames 7 and 8 while the original vortex pair 
shrinks even more. In frame 9 the original vortex pair is no longer 
evident. In frame 10 the two new vortex pairs are fully developed. In 
frames 11 and 12 the larger of the two new vortex pairs rocks to the left. 
Also in frames 11 and 12, the new vortex pair on the left begins to shrink 
toward the concave wall. The larger of the two new vortex pairs is now 
displaced to the left. 

Figure 216 shows a collapsing type of disappearance event. Here, the 
second vortex pair from the right will shrink until it disappears into the 
flow near the concave wall. Frames 1 through 4 show the smallest vortex 
pair present growing to full size. In frames 5 and 6, it twists slightly 
counterclockwise and begins to shrink. Frames 7 through 10 then show the 
rapid collapse into the flow near the concave wall by the small vortex 
pair and then the translation to the left, as it fills the void left 
behind, by the large vortex pair at the right side of the frames. Frames 
11 and 12 also show the center vortex pair enlarging to occupy a larger 
space. 

Figure 217 shows a time history of appearance and disappearance events 
observed over an 80 second interval at De=100 and 115° from the start of 
curvature as observed on the right side of the channel. Figure 218 gives 
magnified views of three regions of Figure 217. The scale on these three 
time lines allows determination of the time intervals between discrete 
events. The time lines show that most appearance events are followed by 


disappearance events, and vice-versa. 
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Frequencies of appearance and disappearance events are given in Tables 
2, 3 and 4 in dimensional and non-dimensional forms. The tabular data 
from Tables 3 and 4 is also shown in Figures 219 to 232. Figures 219, 221 
and 223 present frequencies of appearance and disappearance events versus 
De for angles of 85°, 95° and 105° from the start of curvature 
respectively. Figures 220, 222 and 224 present the same information in 
non-dimensional form as in Figures 219, 221 and 223 respectively. Figures 
225 and 226 present frequencies of appearance and disappearance events 
versus 6 at De=75 and De=100 respectively. Figure 227 shows frequencies 
of total appearance and total disappearance events versus De for angles 
from 85° to 135° after the start of curvature. Figures 228 and 229 
provide curve identification for Figure 227. Figures 230, 231 and 232 
present the same information in non-dimensional form as in Figures 227, 
228 and 229 respectively. Frequencies of appearance and disappearance 
events generally increase with increasing De for Figures 219 through 224 
and 227 through 232. Figure 225 shows that frequencies of appearance and 
disappearance events generally increase with increasing 0 at De=75. 
Figure 226 shows that frequencies of appearance and disappearance events 
generally decrease with increasing 0 at De-100. In all cases, growth type 
appearance events and collapse type disappearance events occur most often. 
In Figures 219 and 225, local minima of total appearance and total 
disappearance event frequencies occur at De-125 and 85? after the start of 


curvature and De-75 and 115? after the start of curvature respectively. 
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V. SUMMARY AND CONCLUSIONS 


Instabilities in curved channel flow were studied over conditions 
ranging from fully laminar to fully turbulent. Primary instabilities in 
the form of regularly spaced Dean vortices develop in the flow due to 
centrifugal effects. After the primary instabilities, secondary 
instabilities develop which affect the behavior of the vortices. These 
secondary instabilities take the form of undulation and twisting of the 
vortices, plus phenomena where vortex pairs appear and disappear. Two 
curved channels with 40 to 1 aspect ratios were used to obtain information 
on these instabilities. 

Time-averaged total pressure, mean velocity, velocity perturbation and 
vorticity component surveys were obtained using a miniature ۴16-06 
pressure probe at De from 50.4 to 424.8 at a streamwise location in the 
second straight section located downstream of the curved portion of the 
channel. At Des75, results show that flow is spanwise uniform. As De 
increases, results provide evidence of vortex pairs located at spanwise 
locations across the measurement plane. At De higher than about 300, flow 
again becomes spanwise uniform. 

Time-averaged and spatially resolved flow correlations were obtained 
uSing a 12 element hot-film probe located 120° from the start of curvature 
at y/d=1.0 (flush with the convex wall) for De from 49.2 to 400.3. Six 
probe elements were used in two different probe configurations with 
element spacings of 6.35 mm and 3.18 mm respectively. In all cases, 
correlation magnitudes decrease as the spanwise distance from the 
reference probe increases. Sign variations of correlations versus z/d 
evidence regions of positive and negative correlation between different 
parts of vortex pair structures. For the first configuration, average 
(over all gages employed) normalized rms voltage increases as De increases 


from De=150 to De=200. Scaled, normalized local rms voltages versus z/d 
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show spatially varying structures at De from 150 to 350. Length scales 
indicate spatially varying structures at De from 100 to 250. 
Qualitatively, the second probe configuration is in close agreement with 
the first configuration. 

Power spectra, streamwise velocity and longitudinal turbulence 
intensity surveys were obtained using a hot-wire probe at De from 50 to 
450. Each survey was obtained at each of five positions in the 
spanwise/radial plane 120° from the start of curvature. Digital and 
analog methods were used to determine streamwise velocity and longitudinal 
turbulence intensity. Results from both methods are in close agreement. 
Streamwise velocity generally increases linearly over the range of De. 
The behavior of velocity for positions near the concave wall varies 
because spatial locations of upwash regions with respect to the concave 
wall vary with De. Longitudinal turbulence intensity increases with De 
higher than 160. Positions near the concave wall have higher intensities 
than the positions at the radial centerline. Observed local frequency 
maxima from low frequency power spectra evidence frequencies of 1 Hz to 2 
Hz, 5 Hz to 6 Hz, 8 Hz, and 20 Hz to 21 Hz for all five probe positions. 
Local frequency maxima also exist at 3 Hz for probe positions away from 
upwash regions. For high frequency power spectra, principle peak 
magnitudes increase with increasing De. Highest peak magnitudes occur at 
positions near the base of the upwash region which evidence twisting [Ref. 
18. 

Frequencies of events where vortex pairs appear and disappear are 
obtained for De from 60 to 150. Time lines showing these events are shown 
at De=100. Still photography time lapse sequences are presented showing 
the collapse of a vortex pair and one type of splitting of a vortex pair. 
These data were obtained for each De from spanwise/radial plane 
observations which were video tape recorded at angles from 85° to 135° 
from the start of curvature. Time line results indicate that should an 


appearance event occur, the likelihood is that a disappearance event will 
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occur next, and vice-versa. Frequencies of total appearance and total 
disappearance events are nearly equal in all cases. Frequencies generally 
increase with increasing De. Growth type of appearance events and 
collapse type of disappearance events occur most often. Highest total 
appearance and total disappearance frequencies observed are at De=125 and 
95º after the start of curvature. Local minima of total appearance and 
total disappearance event frequencies occur at De=125 and 85° after the 


start of curvature and at De=75 and 115° after the start of curvature. 
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APPENDIX A - UNCERTAINTY ANALYSIS 


The uncertainty analysis presented below for sections A. and B. is 
taken from Baun (Ref. 9]. The present study utilizes much of the same 
peripheral equipment and studies flow situations similar to those studied 


by Baun in the same channels. 


A. DEAN NUMBER UNCERTAINTY 
The Dean number is a function of orifice pressure drop. The 
uncertainty for the pressure drop across the orifice is analyzed at 
De=150. 
1. Pressure Drop Uncertainty 
5AP,,.=+0.001 in. H,0 
2. Mass Flow Rate Uncertainty 
6A,,=40.012 in? 
ፅጀ=ቷ0ዐ0 , 003 
ፅሂ=ቷ0 . 02 
6ፀ=ኡ0 . 002 15ጠ/ ፻ቲ* 
ôm=+0.0004 lbm/sec (95% confidence level) 
SA,,=+0.0052 in? 


6De=t4.4 (95% confidence level) 


B. FIVE-HOLE PROBE VELOCITY MEASUREMENT UNCERTAINTY 
Five-hole probe velocity component uncertainties are derived from 
calibration uncertainties and measurement uncertainties. Values at 
De=100.3, y/d=0.5 and z/d=0.0 are listed below. The probe calibration 
uncertainty is for yaw=-10° and pitch=0.0°. 
1. Validyne Transducer Calibration Uncertainty 
ôPranometer £0-001 in. H,0 


5E=+0.0007 Volts 
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§C=+0.000072 in. H,0/Volt (bias error is disregarded because 
pressure differentials are calculated) 


Probe Calibration Uncertainty 
5E=+0.0004 Volts 

5P,=t0.0006 in. H,O 
6P,=t0.0001 in. H,O 
6P,-£0.0005 in. HO 
6P,-£0.0002 in. H,O 

5ጅ5=፡0. 0003 in. H,O 
5P,=+0.00004 in. H50 
5P,=*0.0007 in. H,O 
5P=+0.00015 in. H,O 

05 ۶ رت 5 

.033 0 كة 

5Cbts=+0.045 

6K,=t0.004 degrees”! 

6Kp=t0.004 degrees”?! 

5K,¿=+0.003 degrees jd 

Curved Channel Measurement Uncertainty 
5E=+0.0007 Volts 

6C-£0.000072 1. 456 

5P,=+0.0002 in. H,O 

6P,=+0.0002 in. H,O 

6P;=t+0.0002 in. H,O 

5P,=+0.0002 in. H,O 
6Pç=t0.0002 in. 820 
6P=+0.0001 in. H,O 
6C5,7£0.134 
6+ 41. 
óyaw-t1.33? 


Spitch=+1.4° 
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5(P,-P,)=+0.0004 in. 0 
6U=+0.063 m/s 
ô6U,=+0.062 m/s 
6. 0 742 m/s 


5U,=+0.030 m/s 


C. HOT-WIRE PROBE VELOCITY MEASUREMENT UNCERTAINTY 


The uncertainty associated with the hot-wire probe velocity is derived 


from probe calibration uncertainty and the uncertainty of voltages 


produced by the probe and bridge network. 


1. 


Probe Calibration Uncertainty 
5B=+0.0216 Volts? s“? m-* 

6E =+0.0005 Volts 

6EOC=+0.0005 Volts 

6EOM=+0.0005 Volts 

6AP=+0.0005 in. H,O 

Curved Channel Measurement Uncertainty 
6Eu=+t0.0005 Volts 


5U =+0.048 m/s 
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APPENDIX B - SOFTWARE DIRECTORY 


Programs written in Hewlett Packard BASIC 4.0 language for the HP 300 
Series computers utilized in these studies are listed. A description of 
what the program does, significant user inputs required and significant 


output is provided. 


A. GENERAL PROGRAMS 


DEAN15M Calculates the Dean number as described in Chapter II. 


user input 


e barometric pressure 
e ambient temperature 
e orifice pressure drop 


output 


expansion coefficient Y 
flow coefficient K 
mass flow rate 

channel Reynolds number 
Dean number 

bulk velocity 


B. FIVE-HOLE PRESSURE PROBE PROGRAMS 


C4 Data file of calibration coefficients: pitch, yaw, Coy’ Cop’ 
Cote Cos? and Cots: 

274 Data file of calibration pressures: pitch, yaw, P,, Po, Pz, P,, 
and Ps. 


FIVEHOLE Acquires pressure measurements at each of the five ports 


and directs probe positioning. 


user input 
e date/time 
e estimated dean number 
e orifice size 
e number of spanwise and radial points 
e spanwise and radial resolutions 
e spanwise and radial starting positions 
e wait time for pressure to stabilize 
e number of transducers 
e number of samples per channel 


35 


e number of standard deviations 
e calibration coefficients 
e ambient pressure 


output 


e FIVxxx data file with: y, z, C_, C_ , P,, P, De 
PTUD+>w 。 , pp py 1 = 

e FIVPxxx data file with y, Zz, P4, P5, ,ری‎ ጅረ። Ps, P, De 

e average De 


PADJUST Corrects spatial resolution described in Chapter III. 


user input 


e FIVPxxx data file 
Y len 
Zlen 


output 


* FIVPxxxA data file with y, Zz, C,, C,, P,, P, De 


y' p' 
ORIENT Orients the five-hole probe in the flow such that P =P, at 


yaw=0° and pitch=0° 


user input 


number of samples per channel 
number of standard deviations 
transducer calibration coefficients 
barometric pressure 


output 


mean velocity 
C 
c 


p 


Po, Pz, Por Ps 
VELOCITY1 Computes the three mean velocity components 


user input 


C4 data file 

P4 data file 

number of yaw and pitch points 

FIVPxxxA data file 

pitch angle offset (4.3%) 

yaw angle offset (1.0?) 

y and z direction downwash velocity corrections (0.011, 0.011) 


output 
e Vxxx data file with y, zZz, Per U, Ug, Un, U, 
EA Plots mean streamwise velocity contours 


user input 


e Vxxx data file 
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e number and spacing of contour levels 
output 
e contour plot with level table 
PTOT1 Plots Pamb minus Ptotal contours 


user input 


e Vxxx data file 
e number and spacing of contour levels 


output 
e contour plot with level table 
CCPF1 Calculates CCPF 


user input 


e Vxxx data file 





output 


e PFxxx data file with y, 2, Uccpf， Ux-Uccpf 


UX UCCPF1 Plots streamwise velocity pertubations 


user input 


e PFxxx data file 
e number and spacing of contour levels 


output 
e contour plot with level table 
VORTICITY Computes the three mean vorticity components 


user input 


e Vxxx data file 


output 


D 


e VVxxx data file with y, 57 Og, O., O, 


VORTICITY1 Plots any desired vorticity contour 


user input 
e VVxxx data file 
e number and spacing of contour levels 
e minimum threshold intensity to plot 
e 


type of vorticity to plot (x,y or z) 


output 


e contour plot with level table 
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3DPLOTIxxx Plots a 3-D array of up to seven of the following type of 
plots. Um plots streamwise velocity. م166‎ plots velocity pertubations. 


Pt plots Pamb minus Ptotal. Vx plots .مه‎ Vy plots .مھ‎ Vz plots v All 


z° 
user inputs are identical to those listed in the previously mentioned 
plotting programs with the additional information needed, i.e., the number 


of plots and the identification of all data files to plot. Output is a 3- 


D plot of all graphs: 


ር. 12 ELEMENT HOT-FILM PROBE PROGRAMS 
12FILM1 Acquires raw and conditioned bridge output voltage for six 


active hot-film probes. 


user input 


e active probe numbers 
e De 


output 


e overall sampling frequency 
e individual channel sampling frequency 
+ FDATAxxx data file with active probes, raw and conditioned voltages 


COV2 Analyzes data acquired by 12FILM1 


user input 


e FDATAxxx data file 


output 


mean raw voltages 

conditioned rms voltages 

normalized rms voltages 

Covariances 

correlations 

correlation versus z/d plots for increasing and decreasing z 


D. HOT-WIRE PROBE PROGRAMS 


PROBE 003 Calibration data for the hot-wire probe containing: T 





cal’ 


AS ohr, n, 5), EOC, Ulim and coefficients AO, Al, A2, and A3 
CNT BUF AD Data file used to configure the HP 6944A multiprogrammer 


supplied with the instrument 
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HOTWIRESPM Computes instantaneous streamwise velocities and velocity 


fluctuations. Stores 5000 instantaneous velocities ina data file. 


user input 


sampling frequency 
barometric pressure 
ambient temperature 

bulk velocity 

bridge output zero offset 
Signal conditioner gain 
PROBE 003 data file 
analog mean voltage 
analog rms voltage 


output 


WDATxxxxx data file containing 5000 Us 
channel Re 

De 

digital average Ue 

digital u' 


HOTWIREVEL Computes normalized velocities both digital and analcg. 


Plots Ug versus time. 


user input 


e WDATxxxxx data file 
e U bulk for the given De 


output 


analog U 

analog u’ 

analog and digital u’ 
analog and digital u’/U 
analog and digital ۴ 
plot of Ug 
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response and power spectra plots using the Hewlett Packard 3562A Dynamic 
Signal Analyzer. 


{} codes are user input soft keys with CRT menus. 


APPENDIX C - DYNAMIC SIGNAL ANALYZER OPERATION 


The following sections list sequential operations to acquire frequency 


permanent hard keys shifted to alpha-numeric code. 


A. 


FREQUENCY RESPONSE 


A power spectrum plot from frequency response is obtained as follows: 


Power up the dynamic signal analyzer 
Connect the input BNC connector 
[MEAS MODE] 

(LINEAR RES) 

[SELECT MEAS] 

(FREQ RESP) 

[MEAS DISPLAY) 

{POWER SPEC1} 

[FREQ] 

{FREQ SPAN} 

۲ [ ] - ۱ ٦( 

(Hz) 


This sets up a frequency response from 0.1 to 25 Hz. 


[ SOURCE ] 
(SOURCE OFF) 
[WINDOW] 

(HANN) 

[AVG] 

{STABLE (MEAN) } 
{NUMBER AVGS} 


[5] 


This sets up 5 stable averages. 


[A] 

[SINGLE] 
[COORD] 

{MAG (LOG) } 
{NEXT} 

{LOG X} 
(RETURN) 
[UNITS] 

(P, SPEC UNITS) 
(vé/Hz (PSD)) 
(RETURN) 

(Hz (SEC)) 
(TRACE TITLE) 
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Keys with [] codes are permanent hard keys. Keys with 


Keys with () codes are 


(D) (E) (A) (N) (BD) (2) (BD) [1][0] [0] 
(ENTER) 


This makes the title DEAN = 100. 


[SCALE] 

{X FIXD SCALE} 

] 0 [ ]۰[ )1[ ] ۶ [ ]2[ ]5[ 

{Y FIXD SCALE] 
[0][-][0][0][0] [0][0] [ዐ] [0] [ዐ][ዐ] [ዐ][11[]11] [0] [0] 
[ STATE/TRACE] 


This checks the setup. 


[ STATE / TRACE ] 
[START] 


Wait until the CRT indicates the measurement is complete. 


Ensure the plotter is ready. 


[PLOT] 


POWER SPECTRUM 
A power spectrum plot is obtained as follows: 


Power up the dynamic signal analyzer 
Connect the input BNC connector 
[MEAS MODE] 

{LINEAR RES} 

[SELECT MEAS] 

{POWER SPEC} 

(CH1 ACTIVE) 

(MEAS DISPLAY] 

{POWER SPEC1} 

[ FREQ] 

{FREQ SPAN} 

9۱۱ ۰۱۱۱۱۱۱0۱۲1۱۱۶۱1] [( 
{kHz} 


This sets up a frequency response from 0.1 to 10000 Hz. 


[ SOURCE ] 
{SOURCE OFF} 
[WINDOW ] 

{HANN} 

[AVG] 

{STABLE (MEAN) } 
{NUMBER AVGS} 


[5] 


This sets up 5 stable averages. 


[A] 
[SINGLE] 

[ COORD ] 
{MAG (LOG) } 
{NEXT} 

{LOG X} 
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{RETURN} 
[UNITS] 

{P, SPEC UNITS} 

(V^/Hz (PSD) } 

{RETURN} 

{Hz (SEC) } 

{TRACE TITLE} 

(D) (E) (A) (N) (BL) (2) (BL) [1] [0] (0] 
(ENTER) 


This makes the title DEAN - 100. 


* [SCALE] 
e (X FIXD SCALE) 

[9][-][1][:; ا‎ ۶۶[ 
{Y FIXD SCALE] 
] 0 [ ]۰ [ ] 01 ]0 [ ]0[ ]0[ ]0 [ ] 0[ ]0( ]0[ ]0[ ]0[ ]1[ ] ۶ [ ]1[ ]0 [ ] ۵ [ 
[STATE/TRACE] 


This checks the setup. 


e [STATE/TRACE] 
e [START] 


Wait until the CRT indicates the measurement is complete. 
Ensure the plotter is ready. 


e [PLOT] 
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APPENDIX C - TABLES 


Table 1. 12 ELEMENT HOT-FILM PROBE LOCATIONS 


PROBE # REGION CHARACTERISTIC 
| د‎ [ass | — POSITIVE VORTICITY 
DOWNWASH 
5 
6 







_5.48 | POSITIVE VORTICITY _ 
[5.73 | powwase 
|5ኅ98 | NEGATIVE vortrcrTy | 
0000ا 202ھ‎ | 
EH 
22 





POSITIVE VORTICITY 
5.73 POSITIVE VORTICITY 
DOWNWASH 


12 ያያ. NEGATIVE VORTICITY 


ጋ 
.23 
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Table 2. VORTEX PAIR APPEARANCE AND DISAPPEARANCE OCCURRENCES 


je eae a APPEARANCES DISAPPEARANCES =: 
oen | o f crou | spiss | orot | colt. | nerse engult | cancel | rotat | rime | 


- 
تما 


دا 
ርጋ‏ 
ርጋ‏ 

ርባ 

-3 


0 
دا 





په په 


Table 3. VORTEX PAIR APPEARANCE AND DISAPPEARANCE FREQUENCIES 


——— mmm 


B - APPEARANCES DISAPPEARANCES 
Hee Le pes À imo tec Lene | re | tg | cre | te _ 


a a oa 
125 | 0.08 | 0.00 | 0.08 | 0.08 | 0.00 | 0.00 | 0.60 | 0.05. 
135 | 0.98 | 0.00 | 0.98 | 0.90 | 0.02 | 0.00 | 0.00 | 0.96 | 
] 0.47 | 0.00 | 0.00 | 0.00 | 0.47_ 
ee | 1.22 | 0.00 | 2.22] 3.13 | 0.00 | 0.00 | 0.03 | 1.16 
105 | 1.49 | 0.02 | 2.51] 1.36 | 0.02 | 0.00 | 0.00 | 1.38 
115 | 0.58 | 0.20 | 0.72 | 0.33 | 0.02 | 0.00 | 0.05 | 0.80_ 
385 | 1.70 | 0.02 | 1.72 | 0.50 | 0.35 | 0.12 | 0.65 | 1.72. 
135 | 1.0 | 0.02 | 3.38| 0.57 | 0.35 | 0.9 | 0.23 | 1.22. 
85 [2.96 0.00 [3.01] 2.30 [0.44 [0.00 | 0.02 | 2.76 
[9s |2.31 | 0.12 | 2.48 | 1.97 | 0.08 [0.05 [0.15 | 2.73 | 
[105 |2.02 0.20 | 2.18 | 0.73 | 0.51 | 0.25 | 0.66 | 2.15 
| iis | 2.03 | 0.15 | 2.18 | 0.60 | 0.00 | 0.43 | 0.02 | 2.05 
| 125 | 2.02 | 0.00 | 2.02 | 0.42 | 0.98 | 0.51 | 0.92 | 1.51. 
| 9s | 1.07 | 0.10 | 1.17 | 0.00 | 0.33 | 0.22 | 0.05 | 1.04. 
| 95 | 3.66 | 0.10 | 3.82 | 1.45 | 3.10] 0.02 0.41 | 3.08 


هده هه له > 
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Table 4. VORTEX PAIR APPEARANCE AND DISAPPEARANCE NON-DIMENSIONAL 





DISAPPEARANCES 


Grou | split | Totat | collapse و‎ | wat _ 








D 
rd 
0 
; 
Q 
لکا‎ 
un 





Ul 


.os7]| <o Foo‏ | 0.00 | ده 

Mais. - 71 

781 9.96 | cove | cora eS 

reo [os [70 To.co [ooo [oe [om 
. 094 ee 

عوسی چھا ہہ اود 


125 
135 


ل- 


> 
مم‎ 2 - 
O 


0 
ሀ 


036 


ps 
ps 
Ul 


105 . 106 
105 
. 003 . 107 


.008 .026 
s» | .om | .026 | .012 | .oe1| .101- 
ss [001 | :102 20507077 7 
[osu [o.oo j-o | =o -oae آ تت‎ 
Tas [oas | | خر ره‎ 
65 | -136 | .068 | .142 | .054 | .043 | .039 | .035 | .129. 
205 | .052 | .003 | .055 || .032 | .012 | .014 | .007 | .0451 
150 | 85 | .077 | .005 | .082 | .033 | .024 | .022 | .005 | .064. 
65 | .074 | .007 | .083 | -037 0 


100 T3 





. 107 


E ىا‎ 
O 
ሀ" 
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Table 5. HOT-WIRE PROBE SENSING LOCATIONS RELATIVE TO 


_____ VORTEX PAIR STRUCTURE 
PROBE POSITION lá REGION CHARACTERISTIC 


UPWASH BASE 
ووه سه‎ 


BENEATH VORTEX‏ | 0.10 | 6.60 | 0 پ] 
په 0 وه | اب 
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Table 6. HOT-WIRE SPECTRA SIGNAL CONDITIONER SETTINGS 


SIGNAL CONDITIONING UNITS 


1. Analog E 2. Analog e’ 
Dean Gain Gain 


CE 
Data 


HF plot SEL ۴۱۹۷0۰13 122 362 100 0.1 Hz 1.0 kHz 





* Or as high as possible to not saturate the rms unit 
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Figure 1. Transparent Curved Channel 


49 





Figure 2. Blower Assembly 
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Figure 3. 
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Figure 5. Five-hole Pressure Probe Tip 


58 


x] 
ም Ql ww Oz 00 سا‎ ٢ 


"ëtt? 


٢ وت‎ ፡ 
5 13ባህ፲በህ 170 sal 
ን ۱ 8 ۱ 33n d 
uiu gc —* S 
MIIA 1 
uu DCL p 
| JA Sb X 
—— 4 
با وگ‎ | 
ለ 
MIA 3015 


Suaquinu 0 90 


aqn] 3 —x 


۳4 


p 


Five-hole Pressure Probe Geometry 


Figure 6. 


54 


2- CABELLED 
4^ MALE CONNECTIONS 
.48 cm. O. D. 


95 Cm 





Figure 7. Five-hole Pressure Probe Dimensions and Connections 
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12 Element Hot-film Probe 


Figure 9. 
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Figure 10. 
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Figure 12. Radial/Spanwise Plane Flow Visualization 
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Figure 14. Hot-wire Survey Processing and Data Acquisition 
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Figure 15. 118-35622 Dynamic Signal Analyzer Measurements 
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Figure 16. Hot-wire Probe Locations with Respect to Vortex Pair 
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Figure 17. Total Pressure Contours De=50.4 to De=200.3 
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Figure 18. Total Pressure Contours De=175.2 to De=325.5 
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Figure 19. Total Pressure Contours De=275.2 to De=424.8 
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Figure 20. Streamwise Velocity Contours De=50.4 to De=200.3 
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Figure 22. Streamwise Velocity Contours De=275.2 to De=424.8 
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Figure 24. Velocity Perturbation Contours De=175.2 to De=325.5 
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Velocity Perturbation Contours De=275.2 to De=424.8 
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Figure 26. Streamwise Vorticity Contours De=50.4 to De=200.3 
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Figure 27. Streamwise Vorticity Contours De=175.2 to De=325.5 
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Figure 28. Streamwise Vorticity Contours De=275.2 to De=424.8 
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Figure 29. Radial Vorticity Contours De-50.4 to De-200.3 
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Figure 30. Radial Vorticity Contours De=175.2 to De=325.5 
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Figure 31. Radial Vorticity Contours De=275.2 to De=424.8 
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Figure 32. Spanwise Vorticity Contours De=50.4 to De=200.3 
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Figure 33. Spanwise Vorticity Contours De=175.2 to De=325.5 
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Figure 34. 
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Figure 35. Hot-film Probe Locations with Respect to Vortex Pair 


Structure, De=104.2, 120º Curvature 
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Figure 36. Correlation vs. z/d Probes 2,4,6,8,10,12,wrt 2, 5 
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Figure 37. Correlation vs. z/d Probes 2,4,6,8,10,12,wrt 2, De-100 
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Figure 38. Correlation vs. 2/3 Probes 2,4,6,8,10,12,wrt 2, De=149.9 


86 


Z/D 


™ 

GJ 

CJ 

CU 

1 

LJ 

4 
Ds 
CO 
و‎ 

CU 

O) 

ፎጋ 

CU 

CU 

il 

e 

=) 

8 Y 

< ሏ O -— 


۱۱۵1 10 138803 


Figure 39. Correlation vs. z/d Probes 2,4,6,8,10,12,wrt 2, De=200.7 
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Figure 40. Correlation vs. z/d Probes 2,4,6,8,10,12,wrt 2, De=250.4 
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Figure 41. Correlation vs. z/d Probes 2,4,6,8,10,12,wrt 2, De=299.6 
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Figure 42. Correlation vs. 2/0 Probes 2,4,6,8,10,12,wrt 2, De=350.5 
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Figure 43. Correlation vs. 2/0 Probes 2,4,6,8,10,12,wrt 2, De=399.8 


51 


Z/D 


Figure 44. 


Nu 


DE s 


DE - 358.5 


29926 


\ 


ጂም 
O 
in 
J 
1 
ليا‎ | ~N 
CH ۰ 
C9 
© 
a4 
ዘ 
OB 
© | 。 
0 
wr 
سے‎ 


DE = 
1 0 


3-D Correlations vs. 
De=50.25 to De=399.8 


DE = 


ፍም 


سے 


ae 6; 
ብዥ i 
se ዘ 
x E 
x 
€ S 
NOILUJ 13eleO2 


z/d Probes 2,4,6,8,10,12,wrt 2, 


52 


from 


و 
e‏ 
SO‏ 
LM‏ 
1 
ليا 
للع 
N‏ 
= 
N‏ 
N‏ 
0 
LM‏ 
co‏ 
sr ር...‏ 
co‏ 
( 0 
sr‏ 
CU‏ 
OW‏ 
ዘ‏ 
سے 
c‏ 
SS sr‏ 


NOILB 1313309 


Figure 45. Correlation vs. z/d Probes 2,4,6,8,10,12,wrt 12, De=50.25 
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Figure 46. Correlation vs. z/d Probes 2,4,6,8,10,12,wrt 12, De=100 
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Figure 47. Correlation vs. z/d Probes 2,4,6,8,10,12,wrt 12, De=149.9 
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Figure 48. Correlation vs. z/d Probes 2,4,6,8,10,12,wrt 12, De=200.7 
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Figure 49. Correlation vs. z/d Probes 2,4,6,8,10,12,wrt 12, De=250.4 
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Figure 50. Correlation vs. z/d Probes 2,4,6,8,10,12,wrt 12, De=299.6 
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Figure 51. Correlation vs. z/d Probes 2,4,6,8,10,12,wrt 12, De=350.5 
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Figure 52. Correlation vs. z/d Probes 2,4,6,8,10,12,wrt 12, De=399.8 
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Figure 53. 3-D Correlations vs. z/d Probes 2,4,6,8,10,12,wrt 12, from 
De=50.25 to De=399.8 
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Figure 54. Average Vrms/Vmean vs. De Probes 2,4,6,8,10,12 
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56. Length Scales vs. De Probes 2,4, 6,8, 10,2 
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Figure 57. Correlation vs. z/d Probes 1,2,3,4,5,6, wrt 1, De=49 .2 
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Figure 58. Correlation vs. z/d Probes 1,2,3,4,5,6, wrt 1, De=100.1 
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Figure 59. Correlation vs. z/d Probes 1,2,3,4,5,6, wrt 1, De=150.7 
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Figure 60. Correlation vs. z/d Probes 1,2,3,4,5,6, wrt م1‎ 9 
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Figure 61. Correlation vs. z/d Probes 1,2,3,4,5,6, wrt 1, De=249.6 
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Figure 62. Correlation vs. z/d Probes 1,2,3,4,5,6, wrt 1, De=299.2 
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Figure 63. Correlation vs. z/d Probes 1,2,3,4,5,6, wrt 1, De=350.2 
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Figure 64. Correlation vs. z/d Probes 1,2,3,4,5,6, wrt 1, De=400.3 
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Figure 65. 3-D Correlations vs. z/d Probes 1,2,3,4,5,6, wrt 1, from 
De=49.2 to De=400.3 
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Figure 66. Correlation vs. z/d Probes 1,2,3,4,5,6, wrt 6, De=49.2 
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Figure 67. Correlation vs. z/d Probes 1,2,3,4,5,6, wrt 6, De=100.1 
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Figure 68. Correlation vs. z/d Probes 1,2,3,4,5,6, wrt 6, De=150.7 
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Figure 69. Correlation vs. z/d Probes 1,2,3,4,5,6, wrt 6, De=200.9 
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Figure 70. Correlation vs. z/d Probes 1,2,3,4,5,6, wrt 6, De=249.6 
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Figure 71. Correlation vs. z/d Probes 1,2,3,4,5,6, wrt 6, De=299.2 
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Figure 72. Correlation vs. z/d Probes 1,2,3,4,5,6, wrt 6, 06236 8 
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Figure 73. Correlation vs. z/d Probes 1,2,3,4,5,6, wrt 6, De=400.3 
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Figure 74. 3-D Correlations vs. z/d Probes 1,2,3,4,5,6, wrt 6, from 
De-49.2 to De-400.3 
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Figure 75. Average Vrms/Vmean vs. De probes 1,2,3,4,5,6 
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Figure 76. Vrms/Vmean / Vrms400 
De=400.3 
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Figure 77. Length Scales vs. 
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Figure 78. Digital Mean Streamwise Velocity vs. 
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Figure 79. Analog Mean Streamwise Velocity vs. De, Positions 1,2,3,4,5 
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Figure 80. Digital Longitudinal Turbulence Intensity vs. De, 
Positions 1,2,3,4,5 


128 


DE 


0 ۸ 0 ی‎ 
ALA 





)) EEN 


Figure 81. Analog Longitudinal Turbulence Intensity vs. De, 
Positions 1,2,3,4,5 
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Figure 83. Analog Normalized Longitudinal Turbulence Intensity 
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Figure 85. Peak Frequency vs. De, Position 1 


133 

















(2H) AININDIAS 


Position 2 


De, 
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Figure 88. Peak Frequency vs. 
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Figure 228. Vortex Pair Appearance and Disappearance Frequencies vs. De, 
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Figure 230. Vortex Pair Appearance and Disappearance Non-dimensional 
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Figure 231. Vortex Pair Appearance and Disappearance Non-dimensional 
Frequencies (mfd/U) vs. De, 8=85° to 105° 
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Figure 97. Low Frequency Power Spectrum, De=180, Position 1 
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Figure 104. Low Frequency Power Spectrum, De=50, Position 2 
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Figure 106. Low Frequency Power Spectrum, De 
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Figure 107. Low Frequency Power Spectrum, De=160, Position 2 
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Low Frequency Power Spectrum, 


Figure 141. 
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Figure 142. Low Frequency Power Spectrum, 
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Figure 143. Low Frequency Power Spectrum, De 








ር Se لو ور‎ 617 ۲ SI AS Maa As 
de a TOO AA 2 

| ۳۲۲ |) UE | | Wu AME ر۹‎ | | o OT 
۳ eb EIA یت تست‎ ቅ፦ --፦ لا‎ A |] ات‎ eR 0 > 

| DM d | | II | | 
ٹک‎ CI PES eege የ ዚዚ لت‎ | ጠሌ 
| | ۱ ۲ ۱ ۳ از‎ | | IS , وب‎ 017 | | 
[በ6 ከ6.) 7 ۳ ۱ nO 
TEN E በ | | | ۱ nae i | 
ይ ጠሎ! ie | | | 111 78۹1 | | 
.: ህመ SEES A A A GEN UE ۶[ [E EIE DEM ቤሌ 
| || | | | ۶۹۹۷ ai | | Z A 
SA e. رٹ ا ا ا‎ E 
۱ EET | | | | | l 
EE Nd od 5 le. i [. أ‎ | bt ext 
| [1:1 | | i 
1 | | | 
ር ۶ ARS ۱ 1 188 221 
DII | | 
| eae اا‎ Dc |. 

| ٦ | 

Lu !‏ | د 

| 

| 
| 

| 
)سا 


Position 5 


De=170, 


Low Frequency Power Spectrum, 


Figure 144. 


213 





AAAA J 






نے = 


à 
t 


سے cc መ‏ سے oso as‏ چہ 
zm ላ መጭ መ‏ 

۰ 

። 


۰ 

፥ 

۰ 
= 

E 

| 

| 

۱ 


اډا| 


seo 


一 一 


4 
۱ 
1 
i i 
一 = ہے‎ A 


i 


کے کی نے ہے و ہے نے a A‏ 


د 


— ص س ټس سوہ ا سے سے سد سے —(. سه پس پس پس‎ ጮጮ == መጮ መ ፥ መመ መ سی نس پس مو‎ um 
! ۰ 


A 
«n 


i 
t 
مص س ہے‎ 


ده مهم 


۰ 
۰ 
` 
۰ 
U 
um. ic— m - سد نم س‎ 
፥*=>ሠ መመመ መ سه‎ 


مس .-- 


UL— 


ہے 
١‏ 
4 


( 


Ve تج ۲ سه‎ œ- ል ልመ fe meat mn: a. መ= መሪ ” oa ረ  ፎ<- 
i ۰ 
$ ۰ 
nz cm e کے‎ ss Ss Se 
۰ 


A e 
ፒ 
ሀ] 


t 
i 
— ጨሬመጮመ --“ =- سے‎ 


سه سه ہے سے سو ہم 


Pus 


A ہہ‎ A 


۷ 1 


) 


— .— س سر‎ mm] سے‎ ma 
۱ 
۱ 
t 


N=‏ ۲ تت ېمت 

١ |^ c ክርር 

i ۱ | a ret 

ل | 1 ول 

| | | | 

BOT L OPE 

| | | 1 

Eiere 

| | ۱ ን 

| | ፡ 

[p] ٠ ES 

75۹ داع | | | 

-| Te EP 3 089g 

| | | | 

jo bt c] Bon 

l | ۱ ٦ 

A الو‎ | 

| | | e 

|] | É 

| | | 

| | =| GES 

| Dq | 

| | ۱ ۳ ٣ 

ከ1 ው 

l l | 

E | 1 1 

| | | | 

t lo | ፦- = 

| | | | 

٢ M በበበ مه د‎ 
መ... መመ አሬ 


Position 5 


180, 


De 


214 


Low Frequency Power Spectrum, 


Figure 145. 
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Figure 146. Low Frequency Power Spectrum, De 
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Figure 149. Low Frequency Power Spectrum, De 
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Figure 150. Low Frequency Power Spectrum, De 
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Figure 155. High Frequency Power Spectrum, De 
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Figure 156. High Frequency Power Spectrum, De 
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Figure 159. High Frequency Power Spectrum, De 
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Figure 165. High Frequency Power Spectrum, De 
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Figure 166. High Frequency Power Spectrum, De 


| 


OI 2 1 15602 OS ۵ رح به‎ O AXP* ہے‎ 
ELITS OD TOTSTOET 1,፣ TCI 7 A -1 








۱ )0 77ہ‎ 
۱ 
| 
| 
| 
| 
| 2119. 2. 
لے = دږ‎ 
| 
| BON 
l = oly 
۱ 
۱ 
ا|‎ || ١ | | ۱ 
AAA 
١۱1۱1۱1 | | 
በ. !.0 | | 1. 
[1111 ۱ | | 
۱۲ ۱۱ ۱ | | 
۱۱ ] 1 | 
IIIIII | ۱ ۱۱۱۱ ۱ ۱ ۱ ۱ 0 0 سس‎ ١ ۱ ۱۱۱۱٠ ٠1 | 
|| ተዳ TI == See کت‎ = 
IIIIII | ۱ [111 1 1 | ۱ ۱۱۱ ۱ ۱ 1 | IIIIII Û | unt" ` 
8 OOTI 
DANG 1 2335 كت 230 له‎ 





50, Position 2 


238 


Figure 167. High Frequency Power Spectrum, De 
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Figure 169. High Frequency Power Spectrum, De 
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Figure 170. High Frequency Power Spectrum, De 
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Figure 171. High Frequency Power Spectrum, De 
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Figure 178. High Frequency Power Spectrum, De 
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Figure 193. High Frequency Power Spectrum, De 
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Figure 194. High Frequency Power Spectrum, De 
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Figure 195. High Frequency Power Spectrum, De 
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Figure 196. High Frequency Power Spectrum, De 
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Figure 197. High Frequency Power Spectrum, De 
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High Frequency Power Spectrum, 


Figure 198. 
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Figure 199. High Frequency Power Spectrum, De 
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Figure 200. High Frequency Power Spectrum, De 











MOI Rive zH Bory 6 NEES 
۱۱۲۲ ۲ 1 SEHR, eege 
DEI? ۱ ۲ 1 ۱1 1 locor 
۲۱۱۱۱ + 4 اك‎ E -44141 1-4--ዙ---1 ---- 
7 + 71 ከ ١١٢ TT 111.በ A | 
21111] NMI CRIA ITIN .سل .ابا‎ 
88 ۲۱۱۱۱۱ ۱ 1] | AA 
| RUE ٥ 7 
! ۱۱١١١ ١۱١ | Tee Se 
MTSE 111111 ۹ ۱ 
A ک‎ = 
۳ ۲ jl ۱۱۲۲۲۲ ۱ ۱ al aa 
-TEL E ord MODA 
A 7777 DERA ۷ i 
1 | mcr d ጋ ር 0፡18 ! ۱ 
۱ ۱ ۰۰۹۷ iva O: | Bow 
۱ | | | ZU ٣٠٠٦٠٢ Sea 
1 ٦ | 1018181006. le MERE E d t 
111161 | ۷۱ً ۱۱۲ ۲ لک‎ | 
A. جج و ا‎ 
ET | IIIIII ٠ | ا‎ 
O ما‎ | a Io تمس‎ 
ال‎ d | | IIA ا‎ ROSEN Rad 
10111 1. BM S ۶۰ ۹ A A 
109 [| | | اا‎ | | |. ۲ ۱ "AURI [ono 
۱۱۱۱۱۱ ۱ | ۱ ۱ ۱ | ۱۳۲۲ ۱ [1:11 )] | 107ب‎ 
[11111 Ch ۳ EE E el ہے‎ 
۹ھ‎ |! dd d. od IIIIII I | OS | | ۱۱۱۱۵۲ ۱ L 
ات و‎ ۱1 ۰-۱۱۳۱ ۱ ۰ ۰8۱۱۱۱ > ||| | IA es 
MTT MRI |! 0 eet ۶ ۶ 1 ۳, H | 
INS mp IE t... Ji ic للم اا‎ 15.7... 7 
بانانانا‎ lOO ^N C e, Ee 


Position 4 


400, 


207 


Figure 201. High Frequency Power Spectrum, De 
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Figure 202. High Frequency Power Spectrum, De 
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